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Memorizing Pierre van Baal we will shortly review his life and his scientific achievements. 
Starting then with some basics in gauge field topology we mainly will discuss recent efforts in 
determining the topological susceptibility in lattice QCD. 

Pierre van Baal (Naarden, June 9, 1955 - Leiden, December 29, 2013) was a great theoretician 
many members of the lattice field theory community remember very well. He passed away much 
too early. He was not a ‘latticist’ by himself, but strongly interested in what one can learn from 
lattice theory about fundamental aspects of non-Abelian gauge theories like gluon and quark 
confinement. There are several colleagues who were inspired by his ideas and approaches. The 
author of this contribution and some of his coauthors are grateful to belong to them. 

Pierre started his career with the B.Sc. in Physics and Mathematics in Utrecht. Having received 
his M.Sc. in 1980 he continued with the Ph.D. in Theoretical Physics at Utrecht University, where 
his advisor was Gerard ‘t Hooft. After that in 1984 he moved to Stony Brook first as a Research 
Associate and then as a Fellow in the joint Math/Phys Program. From 1987 to 1989 he became 
a Fellow at the CERN Theory Group. In 1989 he was appointed KNAW-Fellow by the Royal 
Academy of Sciences at University of Utrecht, before he became a full professor in Field Theory 
and Particle Physics at Instituut-Lorentz for Theoretical Physics of the University of Leiden in 
1992. There he was not only a very motivated researcher but also an engaged teacher, even with 
projects for school kids. We all liked him as a very nice, modest person and a good friend of many 
of us. Having been attacked by a serious stroke just after returning from LATTICE ’05 in Dublin, 
he found strong forces to recover and to reestablish his ability to talk, to travel and even to give 
lectures. Unfortunately and tragically, his hope and efforts to come back to research work - as he 
liked it so much - failed. * 1 
But now we feel how much we all miss him. 

The 32nd International Symposium on Lattice Field Theory 

23-28 June, 2014 

Columbia University New York, NY 


* Speaker. 

1 Pierre’s statement in his C.V. (see www.lorentz.leidenuniv.nl/research/vanbaal/DECEASED/HOME/cv.html) 
describes the situation: "I had a stroke (bleeding in the head) on the evening of July 31, 2005. As a consequence of this 

I have accepted that since December 1, 2007 I am demoted to 20% and April 1, 2010 to 10% of a professorship. I could 
still teach (in a modified format), but since October 2008 I can not do it anymore. I can give seminars (twice as slow), 
but doing research (something new) is too difficult." 


(c) Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 
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Figure 1: Pierre van Baal, his field theory lectures book [|l5|] and his work collected by G. ‘t Hooft and C. 
Korthals Altes |[i~6|| celebrated at the Pierrefest in June 2013. 


1. Pierre van Baal’s scientific achievements 


Pierre started his remarkable scientific career in theoretical physics by investigating SU (N) 
gauge fields on a torus in particular with twisted boundary conditions [JT]]. Then he came to his 
“thoughts” on Gribov copies [Q]. Searching for instantons from Monte Carlo generated lattice 
gauge fields he participated in inventing over-improved cooling [^]] which is still in use, as we shall 
report later on. Thinking about improving lattice actions was a further lattice related matter of his 
work 0. His thorough study of multi-instanton solutions and Nahm’s transformation [|51 about 
which several papers appeared over the years [jfij] led him together with his PhD student Thomas C. 
Kraan and partly influenced by papers by Lee and Lu [7|] to his probably most interesting invention 
with the strongest citation impact: periodic instantons (calorons) with nontrivial holonomy 
[j|, ||]. In the following years together with several young (and some senior) collaborators a series 
of papers appeared establishing various properties of those calorons which we want to call in the 
following KvBLL calorons ©0 HI- Some reviews on KvBLL calorons and their relevance can 
be found in his and also other’s ta lk s ©0]. Pierre’s recommendable lectures on field theory 
were published in [15]. His work has been nicely collected in [16] (see Fig. [I]). 


2. Topology, instantons, calorons - a 40 years old story 

Let us first recall some basic facts. Classical as well as path-integral determined quantum 
Euclidean Yang-Mills potentials A^(x) = A a ^(x)T a € su(N c ), tr (T a T b ) = \8 ab , the properties of 
which with their field strength tensor G^ v (x) are defined by the action 

S[A] = — 2 J d x tr (G^vMG^vtr)) > (2-1) 

can be classified by a gauge invariant topological charge 

/ I i 

d x Pt{x), pt{x) = tr (Gnv(x)Gn V (x)), G^ v = —EgvpoGpo . (2.2) 

For finite-action fields in a volume V —> °° the topological charge Q, turns out to be integer-valued, 
because it can be expressed in terms of winding numbers or Pontryagin indices w, of continuous 
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mappings of three-dimensional compact manifolds (surrounding possible singularities of the poten¬ 
tials at finite or infinite x ( ) into the subgroup SU (2), describing “homotopy classes” of the mapping 
^SU{2) = S^\ 

Q t [A] = £>, GZ. (2.3) 

1=1 

The functional Q t [A] is invariant under continuous deformations, which for lattice discretized fields 
holds only if certain smoothness conditions are satisfied. From f d 4 x tr (G uv . ±G^ V ) 2 ] > 0 one 
immediately finds the continuum action S[A] to be bounded from below in each topological sector 


o f 2 

S\A] > —^\Qt[A]\. 

o 


(2.4) 


Gauge field topology became a fundamentally interesting topic for QCD studies almost 40 years 
ago, when classical, topologically non-trivial field configurations, called BPST instantons [ p~7| ], 
were found by solving the (anti) self duality equation G^ v = +<5 UV .. The simplest SU (2) solution 
(with S = \Q,\ = 1) in the singular gauge reads (with 't Hooft’s symbols Pagi = Univ for 
v = 1)2,3, = -p^l = ±8 av , ^2=0) 


a (BPST) 

V 


Iaa„(±) 2 P“ {x z)\ 


(x) =R r\dfiv 


(x-z) 2 ((x-z) 2 + p 2 ) 


(2.5) 


depending on eight modular space coordinates (position z, scale-size p, and global group space 
rotation R ). Note that SU(N c ) solutions are obtained by embedding SU (2) solutions. The sin¬ 
gular gauge instanton potential falls off as x~ 3 at large x. Thus, the Yang-Mills path integral 
[DAexp— S[A] can be semiclassically “approximated” by all possible superpositions of (anti) 
instantons sufficiently distant from each other. Evaluating the integral over leading order quantum 
flucuations around the (anti)instanton configurations, the path integral can be reduced to a parti¬ 
tion function in the modular space of the instanton parameters Jl8[]. This idea led from the dilute 


instanton gas model to an infrared regularized statistical mechanics of an instanton liquid [19]. 
Taking into account Nf fermion flavor degrees of freedom iff with identical mass m, the effect 


of topologically non-trivial configurations like instantons enters through the axial anomaly [20] 


^/ 5 (x) = 2mP(x) + 2Nf p t (x) 


( 2 . 6 ) 


with the topological charge density p r (x) according to Eq. ( |2.2| ) and 


N f 


N f 


f\ X ) = Yj v'/(*)/ t yV/(*)> p ( x ) = L ■ (2-7) 

/=! /=! 

By integrating Eq. ( p.6| ) one gets a relation known as Atiyah-Singer index theorem [21] 

Qt[A]=n + -n^ € Z, (2.8) 


where n + (n ) is the number of zero modes of the massless Dirac operator y- u D 1( [A] with positive 
(negative) chirality on the gauge field background A. A combination of the related Ward identities 
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reads as an identity for the topological susceptibility [ 22 ] 


Xt= \{Qf) 


= / cFx (p t (x)p,( 0)) = - 


N f 


4 m /V -, _ . (2m) 2 

-2 dW/) + — 


( 2 N f ) 


f 


(2 Nf) 2 


d A x (P(x)P(O )), 


1 


m l F n + °(K) 


2N 


f 


(2.9) 


i.e. in full QCD it has to vanish linearly with m 2 in the chiral limit. As we shall see below, 


confirming this limit is still a challenge for lattice QCD. Note that Eq. (2.9) holds also on the lattice 
for Ginsparg-Wilson fermions [23, ^4j] (see below). Applying a 1 /N c expansion, where fermion 
loop contributions become fully suppressed (quenched approximation, i.e. Nf = 0), E. Witten 


(on the basis of current algebra theorems [25]) and G. Veneziano (taking the phenomenological 


spectrum into account [26]) have proposed the relation 


1 


xr= v (Qf 


i 


N f =0 


2N ■ 


F 2 [m 2 , 


+ m 2 — 2m 2 K ] ~ (180MeV) 4 


( 2 . 10 ) 


Therefore, the existence of topologically non-trivial contributions to the path integral leads to the 
solution of the so-called Ua( 1) problem explaining that the tj' meson (of the pseudoscalar flavor 
singlet current) is not a Goldstone boson in the chiral limit and why in nature m^ 3> m K . 

At this place it is worth to note, that the instanton liquid model of the QCD ground state 
describes reasonably well phenomena related to chiral symmetry and Ua( 1) symmetry breaking. 
However, without considering (still unkown) long-range correlations it fails to explain confinement. 
For more information see reviews of instanton physics by T. Schafer and E. Shuryak [|27|] as well 


as by D. Diakonov [|28|], who has passed away also too early. 

Let us turn to the case of non-zero temperature T. The analogous semiclassical treatment of 
the Yang-Mills partition function has been formulated in [ |29| | based on Harrington-Shepard (HS) 
caloron solutions, i.e. ^-periodic instanton chains (1 /T = b) [pG|] 


A" ius) (x) = rja^l <9 v log<l>(jc) 


,(±) 


( 2 . 11 ) 


v-, p 2 Tip 2 sinhf^lx — z\) 

with Ct)fr) — 1 — \ - — - ---- 

hz (x-z) 2 + ( x 4 -Z 4 -kb) 2 b\x — z\ cosh(^|T —z|) — cos (^(x 4 — z 4 )) ’ 


omitting a possible global SU (2) rotation. The topological charge of this solution is 

r b 


2 Jo dM i ^ 3 x P f W = ±1 - 


( 2 - 12 ) 


As for BPST instanton solutions it exhibits trivial asymptotic holonomy, i.e. the untraced Polyakov 
loop at spatial infinity becomes an element of the center Z(2) of SU (2), 


Pexp (ij M(x,x 4 ) dx 4 j 


€Z(2). 


(2.13) 


Today we know that the HS caloron is only a special case of the more general and more compli¬ 
cated KvB LL caloron solution already mentioned above. The SU (N c ) KvBLL caloron in general 
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Figure 2: Left: Slice of the action density at X 4 = const, in a logarithmic scale for a single SU(3) KvBLL 
caloron with non-trivial holonomy showing three static monopole constituents well separated from each 
other. Middle (Right): Localisation of the zero mode for antiperiodic (periodic) boundary conditions for 
the fermion field in the X 4 direction. The figures are taken from Pierre van Baal’s talk at JINR, Dubna in 

1999 0. 


allows non-trivial holonomy, i.e. (f 7JN C ). Such a caloron has Q t = 1 but consists of N c frac¬ 
tionally charged monopole constituents which turn into static (with respect to X 4 ) BPS monopoles, 
if the constituents are sufficiently separated from each other. Because of their selfduality these 
monopoles are often also called dyons. The action or topological charge of the latter are fully 
determined by the eigenvalues of the asymptotic holonomy and constitute together the one- 
instanton action .S'j nst = tin 2 /g 2 of the caloron. Fig. || shows a typical example for the SU (3) case. 
It is obvious that such a configuration cannot be represented as a simple SU (2) embedding. In the 
opposite limit, where the constituents are located near to each other, the action and the topological 
charge density of the KvBLL solution looks very similar to that of a HS caloron (or BPST instan- 
ton), i.e. concentrated within one lump of action and topological charge. Two further properties are 
characteristic for KvBLL solutions. First, the positions of the monopole or dyon constituents are 
given by those locations, where at least two eigenvalues of the local untraced Polyakov loop, i.e. 
the local holonomy , become degenerate 0. Second, the zero mode of the massless Dirac operator 
in a KvBLL caloron background is localized only around one of the constituents (see the middle 
and right panels of Fig. j/j). On which this happens depends on the boundary condition applied to 
the fermion field in the Euclidean time direction [|T(], p~lj ]. 

These properties can be used as a trigger, when detecting calorons and/or dyons from MC gen¬ 
erated thermal lattice gauge field configurations by cooling, 4D APE smearing or overlap operator 
mode expansions (see [I 

SU (3)). The latter series of lattice investigations has led to a simple view of the topological struc¬ 
ture of thermal Yang-Mills fields in terms of (anti)calorons and (anti)dyons. For T < T c , where the 
spatially averaged Polyakov loop is fluctuating around zero, we see all possible dyon constituents 
with equal statistical weight, as one would expect them in a KvBLL caloron with maximally non¬ 
trivial (asymptotic) holonomy. For T > T c , where the Polyakov loop average tends to SU ( N c ) 
center values and where one might expect caloron configurations with holonomies close to such 
values, topological clusters identifiable with corresponding heavy dyon constituents are found sta¬ 
tistically suppressed. As a consequence, on one hand (anti)calorons - containing necessarily heavy 
and light (anti)dyons in this case - are rare, what explains the decreasing topological susceptibility 
with rising temperature. On the other hand, clusters, which can be interpreted (with the triggers 


L 33, 34, 35] for SU (2) pure gauge theory and more recently pq, 37] for 
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mentioned above) as light dyons are abundant. This observation made for SU (2) [34] as well as for 
SU (3) J371 ] supports (non-Abelian) monopole dominance in the deconfinement phase. 

Coming back to the continuum case, the dissociation of KvBLL calorons into dyon con¬ 
stituents has led to the hope to describe quark confinement at T < T c in terms of a liquid of corre¬ 


lated BPS monopoles or dyons. This would realize a picture of instanton quarks | p8| ] studied many 
years ago in detail within the non-linear 0(3)o model [|3^|. Work in such a direction strongly en¬ 
couraged by Pierre van Baal was done over recent years mainly by three groups [ [40| ] (see also talk 
by E. Shuryak). It is worth mentioning that KvBLL calorons carry not only magnetic monopole 
world lines (seen on the lattice within the maximally Abelian gauge) but also extended center vor¬ 


tex structures [35], which seems to provide a bridge to confinement in terms of center vortices, 
too (see the reviews by J. Greensite [|4l|]). Linally, it is worth mentioning that Pierre’s work has 
influenced also a new systematic development of the semiclassical approach within perturbation 
theory called resurgent trans-series expansions |[f2|] (cf. talk by M. Unsal). 


3. How to measure topology on the lattice 

Evaluating the topological charge Q, and correspondingly the topological susceptibility Xt 
as well as identifying topological excitations on the lattice are old issues but remain important 
challenges until today. There are two ways to address this question related to each other via the axial 
anomaly ( |2.6| ) or the index theorem (2.8). The first one expresses the topological charge directly by 
the lattice gluon held strength G^ v . This is easily done with a plaquette loop representation [fO||, 
but the lattice Q t is not an integer, and the corresponding topological susceptibility requires the 
subtraction of a perturbative tail and a proper renormalization p4[]. In combination with various 
methods of stripping off quantum fluctuations i) by cooling - originally invented in order to extract 
approximate multi-instanton solutions |[45||. ii) by 4D APE [|46||. stout [ |47| ] or HYP smearing [fl8||, 
iii) by (inverse) blocking, smoothing or cycling [|fj], or iv) by the gradient how J5(i| ] one ends up 
with smooth lattice gauge held conhgurations, for which (improved) loop dehnitions provide Q t 
values being very close to integers. All these methods applied with a well-defined resolution scale 
allow to reveal the topological structure of the Monte Carlo generated lattice gauge helds in terms 
of clusters of topological charge. 

Alternatives to determine Q t are given by geometric dehnitions of Q t which rely on the ho- 
motopy properties of the gauge held even on the lattice (with torodial boundary conditions). Such 
prescriptions were invented in the past by M. Liischer [51], P. Woit | [52| ] and A. Phillips and D. 
Stone [|53|]. They all provide integer values by dehnition. But due to lattice artifacts on rough 
lattice conhgurations they may yield different numbers. Only sufficiently smoothed helds will pro¬ 
vide a unique answer. Sufficient conditions for such a smoothness exist and can be expressed in 


terms of upper bounds on the action density Q51Q. 

The second approach to determine Q t employs various fermionic dehnitions. The basic obser¬ 
vation is that any lattice Dirac operator obeying the Ginsparg-Wilson relation |fvf|] 


satishes the index theorem 


[ 


y 5 D + Dy 5 = aDy 5 D (3.1) 

Such Dirac operators have been realized within the 


perfect action approach [55], with Neuberger’s overlap operator [57] as well as with domain wall 
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fermions with an extra dimension []58|]. But even for these constructions of Dirac operators holds 
that the topological charge given by their index is not uniquely defined due to lattice artifacts. 

Similarly to applying the gluonic definitions in combination with some smoothing prescription 
one can use a fermionic filtering method by representing the topological charge density p t in terms 
of a finite set of low-lying modes of the choosen lattice Dirac operator D, 

1 N A 

Pt M = tr y 5 (= Ac* - 1 ) = L ('y - 1 ) Vi (X> 75 Vn W • (3-2) 

n =1 


Indeed, comparing this filtering prescription with APE or stout smearing has shown that the number 
of smearing steps can be optimized to a given number N of low-lying modes such that the same 
local topological structures are seen in terms of scale dependent clustering of topological charge 
| [59| l and its varying fractal dimensionality [ |fiC| l. One should expect that this holds also for other 
smoothing methods including also the Wilson or gradient flow. We shall come back to this point 
below (see Section 4.1). 

Other possibilities to determine Q t are given in terms of the anti-Hermitean Dirac operator 0 
PI ] or by counting the index from the spectral flow of the Hermitian Wilson-Dirac operator [|52|]. 
More recently a new fermionic method has been proposed by representing the topological suscep¬ 
tibility in terms of higher moments of scalar and pseudo-scalar currents and spectral projectors 


|, see Section 4.3. 


4. Selected recent lattice results 

4.1 Cooling versus gradient (Wilson) flow 


The old days’ cooling method used for the search for multi-instanton solutions [ fiq ] and more 

|32|, [T 2 L |36|] solves the lat- 


recently to establish the non-trivial holonomy KvBLL calorons [|68|, |32j. [I_2J, 
tice field equations locally (for a given link variable), replaces the old by new link variable, steps 
through the lattice (while the order is not unique), and, if sufficiently often repeated, ends up at 
more or less stable plateau values for the topological charge and action. The method has much im¬ 
proved as over-improved cooling [{|, [S8| stabilizing (multi)instantons or calorons and therefore 
providing extremely stable plateaus at nearly integer Q t values. With an improved lattice represen¬ 
tation of the field strength tensor G^ v [70] one can nicely check the degree of (anti)selfduality e.g. 
by comparing the topological charge with the action in instanton units. 

Recent examples of typical cooling histories for gluodynamics at T > 0 obtained on lattice 
sizes 16 3 x 4 can be seen from Fig. |j. As we have argued in [f7j]] the stability (decay) of plateaus 
for T <T C (T > T c ) can be traced back to the KvBLL caloron structure and to its non-trivial 
(trivial) asymptotic holonomy, i.e. to the dyon constituent mass symmetry (asymmetry). For T > T c 
(anti)selfdual plateaus were occuring very rarely. Searching for with respect to \Q t \ first stable 
plateaus the topological susceptibility Xt for gluodynamics as well as for full QCD with Nf = 2 
clover-improved Wilson fermions is obtained versus T/T c as shown in Fig. f| (for details cf. [7TJ] ). 
It is obvious that Xt behaves much smoother through the crossover of full QCD than passing the 
first order transition in gluodynamics. The behavior of the topological susceptibility at and beyond 
the transition is essential for understanding the mechanism of La( 1) restoration (and in the two- 
flavor case for determining the universality class of the transition in the limit m K —>• 0). It seems to 
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Figure 3: Cooling histories for pure gluodynamics at non-zero temperature (from [[71]]). We show the action 
in instanton units S / Sj nst (blue full lines) and the topological charge Q, (red dashed lines), both represented 
with an lattice-improved field strength tensor, see the text. Left: for confinement at T = 0.887). Right: for 
deconfinement at T — 1.127). 
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Figure 4: Topological susceptibility Xt vs. 777) obtained with over-improved cooling |f7T|| . Left: for 
gluodynamics with lattice sizes 16 3 *4 (red circles), 24 3 *4 (blue up triangles) and 24 3 *6 (green down 
triangles). Right: for full QCD with Nf = 2 clover-improved Wilson fermions for lattice sizes 16 3 * 8 (red 
up triangles) and 24 3 * 8 (blue down triangles). The pion mass is of 0(1 GeV). 
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be too early to draw any final conclusion before taking the chiral limit in a proper way. Therefore, 
it may not wonder, that different groups having recently discussed the issue of Ua (1) restoration 
still come to different conclusions []72|]. 

Compared to (improved) cooling as discussed before, the gradient flow stands theoretically 
on a more sound basis. Proposed and thoroughly investigated by M. Liischer since 2009 [ ^()| l and 
studied also with respect to perturbation theory {73]] it provides an easy controllable manner to 


remove UV fluctuations (cf. his plenary talks at LATTICE 2010 and 2013 |J74[]). It’s flow time 
evolution describing a diffusion process at scale X s — VSt, t = a 2 1 and continuously minimizing 
the action is uniquely defined for an arbitrary lattice field {Up (x) } by solving 


v 4 (*, t) = -g 2 0 [c^S(P (t))] V„ (*, t), Vp (x,0) = Up (x). 


(4.1) 


The physical scale to stop the flow can be efficiently fixed by demanding e.g. 

r(^tr Gp V Gp V )\ t= t, = Tj , i = 0,1 with T 0 = 0.3, T\ = 


(4.2) 
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Simple renormalization properties, in particular in the fermionic sector, and the emergence of topo¬ 
logical sectors at sufficient large diffusion scale are clear advantages of the method. 

However, this does not mean that the previously mentioned cooling or smearing methods have 
to be abandoned. A comparison and mutual optimization of the fermionic filtering method with 
those pure gauge field methods have demonstrated a correspondence between them [|59|]. Indeed, 
the gradient flow can be mapped to cooling in the ensemble average, as C. Bonati and M. D’Elia 
have recently shown J75[]. In the pure gluodynamic case with the standard Wilson plaquette action, 
for given numbers of cooling sweeps n c they have determined the Wilson flow time T, which - in 
the average - provides the same plaquette action. They have estimated the dimensionless flow time 
x on a perturbative ground and found it well satisfied via numerical simulations to be r = n c /3 as 
can be seen from the left panel of Fig. [5j On the right panel of the same figure one nicely sees the 




Figure 5: Left: Flow time T vs. number of cooling steps n c . Right: topological susceptibility Xt vs. 
diffusion scale X s for cooling and gradient (Wilson) flow at various lattice scales. Both figures taken from 



topological susceptibility for cooling and for the Wilson flow to agree completely. Moreover, the 
strong dependence on the lattice spacing at fixed diffusion scale X s becomes obvious. Additionally, 
the authors convinced themselves that cooling and Wilson flow reveal the same local topological 
structure with high confidence. Let us add that recently also Wilson loops have been computed 
with smearing and gradient flow. They were found to agree to a high degree (see [ |76| ] and talk by 
M. Okawa). 


4.2 Exploring the mass dependence of Xt in full QCD 


Only over the last years the expected chiral behavior of the topological susceptibility Xt 


Fjt m n ~ m q{qq) has found a real confirmation from lattice full QCD. Let us briefly sketch some 
recent work in this direction. 


The SINP Kolkata group |78[ ] has employed the standard Wilson gauge and fermion action 
(Nf = 2) at m n > 300 MeV. The topological charge Q t was measured with the blocking-inverse 
blocking (smoothing) method j49[]. An improved ansatz for the density p t was taken. The topo¬ 
logical correlation function for varying volume and quark mass has been studied. A non-negligible 
lattice spacing effect could also be made transparent. From the left panel of Fig. || we observe a 
clear descent towards vanishing Xt with a pion mass squared tending to zero, while on the right 
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L24T48, |i = 5.6 



Figure 6: Left: Xt versus nr K for several volumes. For comparison older data from SESAM-T^L collabo¬ 
ration [f79|] is shown (figure taken from [0D- Right: p, correlation function C t (r) at different quark mass 
values in terms of the Wilson hopping parameter K (from [|78|]). 


panel the space correlation function C,(r) of the topological density shows the expected change in 
sign and a behavior becoming the steeper the smaller the quark mass is. 

A brand-new topology Wilson flow analysis of the ALPHA collaboration has been presented 
at this conference by M. Bruno j8G|] . Since the Wilson flow can be stopped at a well-defined scale, 
the results should be under better control than the previously mentioned ones. The ALPHA col¬ 
laborators studied Nf = 2 lattice QCD with O (a)-improved Wilson fermions and standard Wilson 
gauge action. They investigated the Wilson flow on CLS ensembles with three lattice spacings, 
for m K 6 [190,630] MeV and a lattice extent Lm K > 4. They employed periodic as well as open 
boundary conditions. For the author it came somewhat as a surprise that Q t autocorrelations were 
observed to become weaker with decreasing pion mass 1 . An overall fit to Xt with a ^PT ansatz 
like t\x, = c t\m 2 K + b jy describes the mass dependence sufficiently well (with t\ denoting the 
flow scale according to Eq. (L2)). However, lattice artifacts turned out to be strong and a proper 
chiral limit only possible after the continuum extrapolation is taken (see Fig. [7j). Compared to the 
quenched case also shown in the figure the full QCD result turns out to be quite strongly suppressed 
over the whole range of pion masses studied. 
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Figure 7: Wilson flow estimated topological suceptibility Xt\t=t 0 !2 versus m\ - both in units of t\ - reported 
by the ALPHA collaboration [jgo||. 

Preliminary results of a gradient flow analysis for the topological susceptibility were reported 
'The author thanks S. Mondal for the information that such an observation has been reported also in pT|]. 
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by QCDSF (R. Horsley, G. Schierholz et al) for Nf = 2 + 1 QCD with a tree-level Symanzik 
improved gauge action and (stout smeared) clover-improved Wilson fermions. In this investigation 
QCDSF follows two chiral limit strategies: i) m u = m r j = m s —>• 0 and ii) m u = m r j -y 0, while 
m u + m r j + m s = m = const, with m tuned to its physical value. The corresponding dependence of 
the topological susceptibility on the pion mass can be jointly fitted on the basis of the flavor-singlet 
and flavor-octet Gell-Mann-Oakes-Renner relations (see Fig. |8j). 

Let us conclude this discussion with the comment that the continuum limit for gradient flow 
estimates can be further improved (cf. ta lk s by A. Ramos, S. Sint and D. Nogradi [82, p3t]). 


SU(3) symmetric line m = constant line 




Figure 8: Topological susceptibility fixed at gradient flow scale to versus mr n for two chiral limit strategies 
as explained in the text (communicated by QCDSF, thanks to G. Schierholz). 


4.3 Spectral projector method applied to twisted mass fermions 

Ten years ago, extending an analysis presented in Ref. [|24|l M. Luscher succeeded to pro¬ 


pose a fermionic representation for the topological susceptibility % r in terms of singularity-free 
density chain correlators, which has not to be renormalized []63[] . Treated with spectral projectors 
P m allowing to project onto the subspace of D D eigenmodes below certain threshold M 2 and ap¬ 
proximating them by rational functions R,v/ (see [|64|l) a first computation in pure gluodynamics 
became possible two years later [^5]]. For the valence quarks they used two-flavor clover-improved 
Wilson fermions. The numerical result for ^ quen turned out to be in good agreement with the cor¬ 
responding result [ |84| ] from the index theorem studied with Neuberger’s overlap operator [57] and 
also with the phenomenological value (see Eq. (2.10). To my knowledge, for the first time this 
approach has been applied to compute %, (and the chiral condensate) in full QCD by the ETM 
collaboration [|6^, [67]]. The authors used dynamical Wilson twisted mass fermions with Nf = 2 as 
well as Nf = 2 + 1 +1 flavor degrees of freedom (cf. talks by E. Garcia Ramos and K. Cichy). The 
topological susceptibility has been represented and approximated as 


1 

Xt = 77 


(Tr{R^|) 


F (Tr{y 5 R^y 5 R^}) 


1 


(Tr{y 5 Ri}Tr{y 5 Ri}} 


z 2 

^s 


= v% ( Tr { 75R ^ Tr { 75R w}> = ^2 

where Z(2) random estimators have been used to estimate 

N 


<^ 2 >- 


N 


N 


1 1 

= ^ L ( R v/7VR.v/fl^ R M 75 R Mflr) and = - £ (Ra/Ag 7^ R .v/fl/J • 


k= 1 


N 


(4.3) 


(4.4) 


jt=t 
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Note that the renormalization constants satisfy 0 = 1 and = Q t € Z for N —>• °o, if D is 
a Ginsparg-Wilson operator (e.g. the overlap operator), i.e. 'ff plays conditionally the role of the 
topological charge. For the Wilson twisted mass discretization one can rely on automatic 0(a) 
improvement J85|]. Since the authors used also an improved gauge action they could hope for a 
weak a-dependence of the topological susceptibility. The renormalization constants Z 5 ,Zp can be 
taken also from other ETMC evaluations [[8h||. The “topological” charge '(> turned out nicely 
Gaussian-like distributed. In Fig. |9] selected results are shown. The left panel demonstrates the 
projector method computation of the Z-factor ratio to be consistent with that of Refs. [|S6||. The 
right panel shows the topological susceptibility as a function of the quark mass. Within the error 
bars, which are still quite large, the different lattice scale results more or less agree, indeed. In any 
case the behavior of Xt is seen to be compatible with a linear decrease with the quark mass towards 
the chiral limit. The slope of this curve allows also to determine the chiral condensate. The result 
of this estimate was consistent with that of other methods. Finally, in the quenched limit ^ r quen 
came out in good agreement with Eq. ( 2.10 ) (cf. E. Garcia Ramos’ talk). 
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Figure 9: Left: The ratio of renormalization constants 0 versus threshold mass Mr for Nf = 2. In the right 
most part of the panel results of alternative computations are shown [^6|]. Right: The Nf = 2+1 + 1 result 
for Xt in units of the Sommer scale ro | |87| | versus renormalized quark mass )Ir for three lattice spacings. 
Figures are taken from [j57||. 


4.4 Comparing various methods to determine Q t and Xt 

The ETM collaboration has made a joint effort J88[ ] to compare various methods to compute 
the topological susceptibility within the framework of Nf = 2 twisted mass fermions and tree-level 
Symanzik improved gauge action. All computations were done on the same set of configurations 
for a pion mass value m n = 300 MeV, with a lattice spacing a = .081 fm and linear lattice size L = 
1.3 fm. Without taking the continuum limit one should not really expect a full agreement even in 
case the methods are really equivalent and correctly established. The authors compared fermionic 
definitions of the topological charge (index of the overlap operator, Wilson-Dirac operator spectral 
flow (SF), spectral projector method (SP)) with gluonic field strength (FT) definitions applying 
various versions and stages of gradient flow (GF), cooling, and APE/HYP smearing. In Fig. [l0| 
we see how the different methods to determine Q, are correlated. On the right hand side the 
colored correlation scale from bottom (blue = weak correlation) to top (red=strong correlation) 
is given. Most of the versions are strongly correlated among each other, except the gluonic FT one 
without any removal of perturbative fluctuations. This is a well-known fact since the early days 
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of topological lattice investigations. That also the spectral projector results are weakly correlated 
to the other ones does not come as a surprise, too, because of the stochastic determination of the 
‘charge’ c tf in accordance with Eq. (£Oj). A comparison of the resulting ft values shows that all 
methods provide results in the same ball park. This holds also for the spectral projector method, 
provided the spectral threshold Mr is taken high enough. For details we refer to [88]. The outcome 


naive FT basic cooling 30 
naive FT basic cooling 10 
impr. FT basic cooling 30 
impr. FT basic cooling 10 
naive FT impr. cooling 30 
naive FT impr. cooling 10 
impr. FT impr. cooling 30 
impr. FT impr. cooling 10 
naive FT APE30 
naive FT APE10 
impr. FT APE30 
impr. FT APE10 
w impr. FT HYP30 

| impr. FTHYP10 

impr. FT nonSmear 
£ GF flow time 3t0 

t-j GF flow time 2t0 

GF flow time to 
spec. proj. M 2 =0.0015 
spec. proj. M 2 =0.0010 
spec. proi. M=0.0004 
spec. proj. M=0.0000355 
SF HYP5 s=0.0 
SF HYP5 s=0.5 
SF HYP1 s=0.0 
SF HYP1 s=0.75 
index HYP1 s=0 
index nonSmear s=0 
index nonSmear s=0.4 
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Figure 10: Correlation matrix between different methods to determine the topological charge (from |^8[j). 


of the ETMC comparison certainly allows the conclusion that the well-controllable gradient flow 
method at the first place, but also cooling or smearing methods, if adapted with analogous criteria to 
fix the diffusion scale, are optimal (also from the computational point of view) in order to determine 
topological properties of lattice QCD. 


5. Conclusions 


First of all I have to apologize for not having discussed many issues, which have been touched 
during this symposium and would have been of interest here - as they would for Pierre van Baal. 
The following topics belong to the problems within the range of this short review but could not be 
covered here: 

- the recent rj' — rj mixing results [f*Tj|, which became possible due to various powerful noise 
reduction techniques [fj()||. 


the £4(1) symmetry restoration puzzle mentioned already above [72], 


- the use of open boundary conditions suppressing HMC’s autocorrelation for Q t [91, 80] 

(cf. talk by G. Me Glynn), 

- the simulation of 0-vacua with Langevin techniques or dual variables 
(talks by L. Bongiovanni, T. Kloiber), 

■ considerations with fixed topology (talks by A. Dromard, H. Fukaya, U. Gerber, J. Verbaarschot), 

■ ongoing discussions about the vacuum structure and topological excitations (talks by N. Cundy, 
P. de Forcrand, M. Hasegawa, M. Ogilvie, A. Shibata, H.B. Thacker, D. Trewartha, M. Unsal), 


■ phase structure at differing masses ([92] and talk by S. Aoki), 

■ topology in related theories as G 2 Yang-Mills theory and N = 1 SUSY on the lattice 
(|p3|| and talk by P. Giudice), 


- effects in QCD caused by magnetic background fields [94]. 
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Let us briefly summarize. Investigations of KvBLL caloron and dyon gas models with non¬ 
trivial holonomy as initiated by Picric van Baal are still interesting and encouraging for better 
describing (de)confinement within the finite temperature setting of QCD. Even more, they may 
pave a way to improve systematically the semiclassical approach. The computation of the topolog¬ 
ical susceptibility with new methods (gradient flow, spectral projector method) is on a promising 
way. In any case one has to keep track of lattice artifacts and to study the continuum limit. 
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